home *** CD-ROM | disk | FTP | other *** search
/ 64'er / 64ER_CD.iso / s85xx / s8507b.d64 / runge kutta (.txt) < prev    next >
Commodore BASIC  |  1995-03-30  |  507b  |  23 lines
  1. 100 PRINT"[147] RUNGE-KUTTA-VERFAHREN OHNE LOK. FEHLER [146]"
  2. 110 PRINT"DGL AM PROGRAMMENDE EINGEBEN!"
  3. 120 INPUT"SCHRITTWEITE  ";H
  4. 130 INPUT"ANFANGSWERT   ";Y0
  5. 140 INPUT"INTERV.GRENZEN";A,B
  6. 150 DEFFNR(X)=INT(X*1E4+.5)/1E4
  7. 160 PRINT" X"," Y"
  8. 170 YK=Y0
  9. 180 XK=A
  10. 190 GOSUB250
  11. 200 YL=YK+(H/6)*(L1+2*L2+2*L3+L4)
  12. 210 PRINTFNR(XK),YK
  13. 220 XK=XK+H:YK=YL
  14. 230 IFXK<=B+HTHEN190
  15. 240 END
  16. 250 X=XK:Y=YK:GOSUB300:L1=F
  17. 260 X=XK+H/2:Y=YK+(H/2)*L1:GOSUB300:L2=F
  18. 270 X=XK+H/2:Y=YK+(H/2)*L2:GOSUB300:L3=F
  19. 280 X=XK+H:Y=YK+H*L3:GOSUB300:L4=F
  20. 290 RETURN
  21. 300 F=Y*X-2*X  :REM <=> Y'=... _DGL
  22. 310 RETURN
  23.